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Context. The Large Magellanic Cloud (LMC) is to date the only normal external galaxy that has been detected in high-energy gamma rays. High- 
energy gamma rays trace particle acceleration processes and gamma-ray observations allow the nature and sites of acceleration to be studied. 
Aims. We characterise the distribution and sources of cosmic rays in the LMC from analysis of gamma-ray observations. 
' . . Methods. We analyse 1 1 months of continuous sky- survey observations obtained with the Large Area Telescope aboard the Fermi Gamma-Ray 
^ Space Telescope and compare it to tracers of the interstellar medium and models of the gamma-ray sources in the LMC. 
k> ^^^^1^^- The LMC is detected at 33cr significance. The integrated >100 MeV photon flux of the LMC amounts to (2.6 + 0.2) X 10"^ ph cm"^ s"^ 
^\ which corresponds to an energy flux of (1.6 + 0.1) x 10~^^ erg cm~^ s~\ with additional systematic uncertainties of ^ 16%. The analysis reveals 
I the massive star forming region 30 Doradus as a bright source of gamma-ray emission in the LMC in addition to fainter emission regions found 
in the northern part of the galaxy. The gamma-ray emission from the LMC shows very little correlation with gas density and is rather correlated 
to tracers of massive star forming regions. The close confinement of gamma-ray emission to star forming regions suggests a relatively short GeV 
cosmic-ray proton diff'usion length. 

Conclusions. The close correlation between cosmic-ray density and massive star tracers supports the idea that cosmic rays are accelerated in 
massive star forming regions as a result of the large amounts of kinetic energy that are input by the stellar winds and supernova explosions of 
massive stars into the interstellar medium. 

Key words. Acceleration of particles - cosmic rays - Magellanic Clouds - Gamma rays: observations 



1. Introduction 

Since the early days of high-energy gamma-ray astronomy, it has 
been clear that the gamma-ray flux receive d on Earth is dorn - 
inated by emission from the Galactic disk (IClark etalJ 119681) . 



This emission is believed to arise from cosmic-ray interac- 
tions with the interstellar medium, which at gamma-ray energies 
^100 MeV are dominated by the decay of tt^ produced in col- 
lisions between cosmic-ray nuclei and the interstellar medium 
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(lPollack&Faziolll963h . Further contributions are from cosmic- 
ray electrons undergoing inverse Compton scattering off inter- 
stellar soft photons and Bremsstrahlung losses within the inter- 
stellar medium. Gamma-ray observations thus have the potential 
to map cosmic-ray acceleration sites in our Galaxy, which may 
ultimately help to identify the sources of cosmic-ray accelera- 
tion. 

Nearby galaxies have some advantages as targets for stud- 
ies of cosmic-ray physics and have the advantage of being 
viewed from outside, and so line of sight confusion, which com- 
plicates studies of emission from the Galactic disk, is dimin- 
ished. This advantage is, however, somewhat offset by the lim- 
itations by the angular resolution and sensitivity of the instru- 
ment. The Large Magellanic Cloud (LMC) is thus an excellent 
target for studying the link between cosmic-ray acceleration and 
gamma-ray emission since the galaxy is nearby {D ^ 50 kpc; 
[Mat sunaga et al. 2009; Pietrzynski et al. 2009), has a large an- 
gular extent of ~ 8°, a nd is seen at a small inclination an - 
gle of / ^ 20° - 35° ( Kim et al.1 [T998[: IVan der Marei[l200^ . 
which avoids source confusion. In addition, the LMC is rela- 
tively active, housing many supernova remnants, bubbles and 
superbubbles, and massive star-forming regions that are all po- 
tential sites of cosmic -ray acceleration (ICesarsky & Montmerl^ 
ll983l:lBiermann[l2004l:lBinns et al.ll2007L 

The EGRET telescope aboard the Compton Gamma-Ray 
Observatory (CGRO, 1991-2000) was the first to detect the 
LMC (Sree kumar et al.lll9 92), which has remained the only nor- 
mal galaxy besides our Milky Way that has been seen in high- 
energy gamma rays. Due to EGRET 's limited angular resolu- 
tion and limited sensitivity, details of the spatial structure of the 
gamma-ray emission could not be resolved, yet the observations 
showed some evidence that the spatial distribution is consistent 
with the morphology of radio emission. The observations also al- 
lowed determining the integral gamma-ray flux from the LMC, 
which was shown to be consistent with expectations based on 
a mo del using the princi ples of dynamic balance and contain- 
ment ([Fichtel et al.[[l99l[) . From the agreement with this model, 
[Sreekumar et al.[ ([1992 *) conclude that the level of cosmic rays in 
the LMC is comparable to that in our Galaxy. 

The Large Area Telescope (LAT) aboard Fermi is now for the 
first time providing the capabilities to go well beyon d the study 
of the integrated gamma-ray flux from the LMC CDigel et al.[ 
[2000[: [^idenspointner et al.[[2007[) . Thanks to its excellent sen- 
sitivity and good angular resolution, we are able to spatially 
resolve the structure of the gamma-ray emission, allowing us 
to provide a detailed mapping of the cosmic-ray density in the 
galaxy. In the spectral domain, our data allow for detailed com- 
parisons with cosmic-ray interaction models and can be used 
to search for spectral variations over the galaxy. In the tim- 
ing domain, we are able to assess contributions from flaring 
sources and can search for pulsations from energetic p ulsars in 
the galaxy that may c ontribute to the total emission ([Harding[ 
ll98ll:lHartmann et al.|[T993[) . 

In this paper we present our first in-depth analysis of the 
LMC galaxy based on 1 1 months of continuous sky survey ob- 
servations performed with FermifLAI. We put particular em- 
phasis on determining the spatial distribution of the gamma-ray 
emission, which, as we will show, reveals the distribution of cos- 
mic rays in the galaxy. We also search pulsations from energetic 
pulsars that may contribute to a non-negligible level to the over- 
all emission. 



2. Observations 

2.1. Data preparation 

The characteristics and p erformance of t he LA T aboard Fermi 
are described in detail by lAtwood et al.[ ([2009 ^). The data used 
in this work amount to 274.3 days of continuous sky survey ob- 
servations over the period August 8 2008 - July 9 2009 dur- 
ing which a total exposure of ~ 2.5 x 10^^ cm^ s (at 1 GeV) 
has been obtained for the LMC. Events satisfyin g the standard 
low-background event selection ("Diffuse" events; Atwood et al.[ 
2009) and coming from zenith angles < 105° (to greatly re- 
duce the contribution by Earth albedo gamma rays) were used. 
To further reduce the effect of Earth albedo backgrounds, the 
time intervals when the Earth was appreciably within the field 
of view (specifically, when the centre of the field of view was 
more than 47° from the zenith) were excluded from this analy- 
sis. Furthermore, time intervals when the spacecraft was within 
the South Atlantic Anomaly were also excluded. We further re- 
stricted the analysis to photon energies above 200 MeV; below 
this energy, the effective area in the "Diffuse class" is relatively 
small and strongly dependent on energy. 

For the analysis we selected all events within a rectan- 
gular region-of-interest (ROI) of size 20° x 20° centred on 
(q^j20oo, ^J20oo) = (05M7"^36^ -69°0r48"f| and aligned in equa- 
torial coordinates. All analysis was performed using the LAT 
Science Tools package, which is available from the Fermi 
Science Support Center, using P6_V3 post-launch instrument re- 
sponse functions (IRFs). These take into account pile-up and ac- 
cidental coincidence effects in the detector subsystems that were 
not considered in the definition of the pre-launch IRFs. 



2.2. Data selection and background modelling 

At the Galactic latitude of the LMC (b ^ -33°), the 
gamma-ray background is a combination of extragalactic and 
Galactic diffuse emissions and some residual instrumental 
background. The extragalactic component comprises resolved 
sources, which of ten can be associated with known blazars 
dAbdo et al.1 [2009b) . and a diffuse component, which is at- 
tributed to unre solved sources and eventually intr insically dif- 
fuse processes ("Strong et al.[ [2004a[: [Dermed [2007 ). The latest 
LAT collaboration-internal source list that has been derived from 
1 1 months of survey data (comparable to the data volume used 
in this analysis) contains 542 sources for latitudes \b\ > 30°, 
corresponding to a source density of 86 sources sr"^ Our ROI 
covers a solid angle of 0.12 sr, so we expect about 10 resolved 
background sources in our field. Among those, 1-2 should spa- 
tially overlap with the LMC if we assume the galaxy's diameter 
is between 8° and 10°. 

Within the ROI but outside the boundaries of the LMC we 
found a total of 6 significant point sources, compatible with 
expectations (cf. Table [T]). For all of them we found counter- 
parts in the CRAT ES catalogue of flat-spectrum ra dio sources 
(iHealev et al][2007 ) using the procedure described in [Abdo et al.[ 
([2009a[) . making them good candidates for being background 
blazars. Searching for background sources within the LMC 
boundaries is more difficult due to possible confusion with local 
emission maxima of the galaxy's diffuse emission. Background 
blazars may however reveal themselves by their flaring activ- 
ity, hence we searched for any excess emission above the time 
averaged level on a monthly basis (cf. section [24l) . We found ev- 
idence for a flaring source during month 4 of our dataset (MJD 



corresponding to (l,b) = (279.73°, -33.54°) 
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Fig. 1. Gaussian kernel (cr = 0.2°) smoothed counts maps of the region of interest (ROI) in a true local projection before (left) 
and after subtraction of the background model (right) for the energy range 200 MeV - 20 GeV and for a pixel size of 0.1° x 0.1°. 
Overlaid is the N(Hi) contour of 1 x 10^^ H cm"^ of the LMC to indicate the extent and shape of the galaxy. The boxes show the 
locations of the 6 point sources that have been included in the background model. The right panel has a true dynamic range from 
-46 to +248 counts deg"^ that has been expanded for display to cover the full dynamic range of the residuals that are shown in 
Fig. HI 



Table 1. Point sources included in the background model. 



Name 




TS 


Q^J2000 




CRATES J050755- 


-610441 


118.6 


05^07^^55^ 


-6r04'43" 


CRATES 105 1643- 


-620706 


269.3 


05M6"^45^ 


-62°07'05" 


CRATES J055842- 


-745904 


53.5 


05^58^^46^ 


-74°59'05" 


CRATES J060106- 


-703606 


52.5 




-70°36'09" 


CRATES J063542- 


-751615 


90.0 


06^35^^46^ 


-75n6'16" 


CRATES J070027- 


-661041 


123.4 


07hOO"^3P 


-66° 10' 45" 



Note to th e table : TS is a measure of the detection significance of the 
source fcf. l23^ . 



54777.8 - 54808.2) near 30 Doradus, and thus we excluded the 
data within this time interval from our analysis. This results in a 
dataset that corresponds to 248.7 days of continuous sky survey 
observations during which a total exposure of ~ 2.3 x 10^^ cm^ s 
(at 1 GeV) was obtained for the LMC. 

We then modelled background gamma-ray emission within 
the ROI using components for the diffuse Galactic and the ex- 
tragalactic and residual instrumental backgrounds and the 6 
blazars. The Galactic component was based on the LAT standard 
diffuse background model gll_iem_v020 for which we kept the 
overall normalisation as a free parameter. The extragalactic and 
residual instrumental backgrounds were combined into a single 
component which has been taken as being isotropic. The spec- 
trum of this component was determined by fitting an isotropic 
component together with a model of the Galactic diffuse emis- 
sion and point sources to the data. Also here we left the overall 



^ The mo del can be d ownloaded from 
|http://fermi.gsfc.nasa.gov/ssc/data/access/lat/Back groundM odels.htmlt 



normalisation of the component as a free parameter. The 6 back- 
ground blazars were modelled as point sources with power-law 
spectral shapes. The positions of the blazars were fixed to those 
given in the CRATES catalogue ( Healey et al. 2007) and are 
given in Table [T] The flux and spectral power-law index of each 
source were left as free parameters of our background model and 
their values were determined from likelihood analysis. 

2.3. Spatial distribution 
2.3.1. Counts map 

To investigate the spatial distribution of gamma-ray emis- 
sion toward the LMC we first binned all photons into a 
counts map of size 20° x 20° centred on (0^12000, ^J20oo) = 
(05^17"^36^ -69°0r48") and aligned in equatorial coordinates. 
Figure[T] shows the counts map before (left panel) and after (right 
panel) subtraction of the background model. 

The background subtracted map confirms that diflTuse 
Galactic and isotropic backgrounds as well as the 6 background 
blazars are properly removed by our treatment. The only re- 
maining feature is extended emission that is spatially confined 
to within the LMC boundaries which we trace by the iso col- 
umn density contour A^h = 1 x 10^^ H cm"^ of neutral hydrogen 
in the LMC (Kim et al.ll2005l) . The total number of excess 200 
MeV - 20 GeV photons above the background in the LMC are43 
amounts to ~ 1550 counts whereas the background in the same 
area amounts to ~ 2440 counts. With these statistics, the ex- 
tended gamma-ray emission from the LMC can be resolved into 
several components. The brightest emission feature is located 

^ We use a square region 5.5° x 5.5° centred on (Qfj2ooo. <^j2ooo) = 
(05^30"^, -68°30') to extract these numbers. 
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Fig. 2. Profiles of the observed number of counts along Right 
Ascension (top) and Declination (bottom) of the LMC region for 
the energy range 200 MeV - 20 GeV. The solid lines show the 
corresponding profiles of the background model. The red arrow 
indicates the location of 30 Dor. 



near (aj2ooo, ^j20oo) ~ (OSMO"", -69° 15'), which is close to the 
massive star-forming region 30 Doradus (30 Dor) that houses the 
two Crab-like pulsa r s PSR J0540-6919 a nd PSR J0537-6910 
(ISeward et al.1 1 19841 : [Marshall et aP 1 19981) . Excess gamma-ray 
emission is also seen toward the north and the west of 30 Dor. 
These bright regions are embedded into a more extended and 
diff'use glow that covers an area of approximately 5° x 5°. To 
further illustrate the emission structure we show profiles of the 
counts map in Fig. [2] along aj2ooo (top panel) and ^j2ooo (bot- 
tom panel) and we provide an adaptively smoothed image of the 
30 Dor region in Fig.O 

2.3.2. Geometrical models 

Out next step was to assess the spatial distribution of the LMC 
emission using simple parametrised geometrical models of the 
gamma-ray intensity distribution. We assumed power-law spec- 
tral distributions for all models and kept the total flux and power 
law index as free parameters. We adjusted the spatial and spec- 
tral parameters of the models using a binned maximum like- 
lihood analysis with spatial bins of 0.1° x 0.1° and 60 loga- 
rithmically spaced energy bins covering the energy range 200 
MeV - 20 GeV. For each model we computed the so-called Test 
Statistic (TS), which is defined as twice the diff'erence between 
the log-likelihood Xi that is obtained by fitting the model on 
top of the background model to the data, and the log-likelihood 
Xo that is obtained by fitting the background model only, i.e. 
TS = 2(Xi - Xo). Under the hypothesis that the background 
model satisfactorily explains our data, TS follows the distri- 
bution with p degrees of freedom, where p is the number of ad- 
ditional free parameters in the maximisation of Xi with respect 
to those used in the maximisation of Xo ( Cash 1979). For exam- 
ple, for a single source on top of the background model with free 
position, flux and spectral power law index we have p = A. 

First, we examined whether the gamma-ray emission from 
the LMC can be explained with a combination of individ- 
ual point sources. For this purpose we added successive point 
sources to our model and optimised their locations, fluxes and 
spectral indices by maximising the likelihood of the model. We 
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Fig. 3. Zoom into a 3° x 3° large region of the background 
subtracted counts map around the central star cluster R 136 
in 30 Dor. The c ounts map has been adaptively smoothed 
dEbeling et al.ll2QQ6 ) with a signal-to-noise ratio of 5 in order to 
reveal significant structures at all possible scales while suppress- 
ing the noise that arises from photon counting statistics. Stars 
show the locations of R 136 and of the pulsars PSR J0540-6919 
and PSR J0537-6910. The circles show the 95% containment 
radius of sources SI, S3, and G2 (see text). For G2 we added the 
Gaussian cr of the source to the error radius. We also show the 
95% containment radius of a flaring source that appeared during 
month 4 of the analysed data set. The diamond shows the loca- 
tion of the possible counterpart RX J0546. 8-6851 of the flaring 
source (see section [24b . 



stopped this procedure once the TS improvement after adding 
a further point source dropped below 25. This happened after 
we added a 6th point source to our model which resulted in a 
TS improvement of only 21.9, corresponding to a detection sig- 
nificance of 3.7(T for this source. Table [2l provides the maximum 
likelihood positions and fluxes for the 5 significant point sources, 
and gives also the TS improvement (ATS) that is achieved in 
each of the successive steps. Adding up ATS for all sources pro- 
vides a total TS of 1089.3 for this model, corresponding to a 
detection significance of 31.5(T (/? = 20). We refer to the point 
source model as PS. 

Second, instead of using point sources we repeated the pro- 
cedure with 2D Gaussian shaped intensity profiles to build a geo- 
metrical model that is more appropriate for extended and diff'use 
emission structures. We again stopped the successive addition 
of 2D Gaussian shaped sources once the TS improvement after 
adding a further source dropped below 25. This occurred after 
two 2D Gaussian shaped sources have been added to the model. 
Table[3]lists the maximum likelihood positions, Gaussian widths 
cr and fluxes for the two significant 2D Gaussian sources. Similar 
to the point sources, the column ATS quotes the TS improve- 
ments in each of the successive steps. The total TS amounts to 
1 122.6 for this model, corresponding to a detection significance 
of 32.7(T {p = 10). We refer to this model from now on as 2DG. 
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Table 2. Parameters of the point source model PS. 



Name 


ATS 


Q^J2000 


<^J2000 


r95 


Flux 


SI 


798.1 


05^38.9"^ 


-69n7' 


A' 


11.8 + 1.8 


S2 


163.8 


05^30.2"^ 


-67n6' 


IT 


7.6 + 2.1 


S3 


66.5 


05^27.7"^ 


-68°58' 




4.0+1.7 


S4 


33.1 


05Ml.2"^ 


-67n5' 


25' 


3.9 + 2.3 


S5 


27.8 


04^55.8"^ 


-69°21' 


\V 


2.6 + 1.4 



Note to the table: Columns are (1) source name, (2) TS increase after in- 
cluding the respective point source in the model, (3) Right Ascension of 
source, (4) Declination of source, (5) 95% containment radius, and (6) 
integrated > 100 MeV photon flux in units of 10~^ ph cm~^ s~^ assum- 
ing a power law spectral slope with free index for each source. Quoted 
uncertainties are statistical only. 

Table 3. Parameters of the 2D Gaussian sources model 2DG. 
Columns are identical to those in Table [21 except for column 6 
which gives the width cr for each 2D Gaussian source. 



Src. ATS Qfj2ooo ^moo r95 cr Flux 

Gl 1000.9 05^26.0"^ -68°16' 20' 73' + 5' 19.6 + 2.2 
G2 121.7 05^38.8"^ -69^18' 7' 12' + 4' 8.5 + 2.2 



Adding a 3rd 2D Gaussian shaped source improved the TS 
by only 21.9, corresponding to a detection significance of 3.5cr 
for this source (p = 5). We note that this 3rd 2D Gaussian 
source is located at (ckj2ooo, ^J20oo) = (05^17. 8"", -72°320 with 
a 95% containment radius of 18', and that it formally is con- 
sistent with a point source (cr = W ± 13'). Using our stan- 
dard source association procedure (lAbdo et al.|[2009a ) we found 
the blazar CRATES J05 1636-723707 as a possible counterpart 
of this source, suggesting the presence of a background blazar 
within the boundaries of the LMC (see also section [231) . 

To illustrate how well the models fit the data we show in 
Fig. m the model counts map and the residual counts that are left 
in the ROI after subtracting the PS or 2DG model from the data. 
The colour scale and dynamic range of the residual maps has 
been chosen identical to that of the right panel of Fig.[T]to allow 
the comparison of the residuals before and after subtraction of 
the LMC model. Apparently, only very few residual counts are 
left in the LMC area after subtraction of either of the two mod- 
els from the data (however a significant negative residual at the 
position of S 1 indicates an overestimation of the flux in the PS 
model in this area). None of the peaks in the residuals is statisti- 
cally significant (ATS < 25) as a point source (for the PS model) 
or a 2D Gaussian shaped extended source (for the 2DG model). 

The PS and 2DG models have in common that they 
both comprise a strong source located at (ckj2ooo, ^J20oo) ~ 
(05^39^, -69° 170 in the 30 Dor region (SI for PS and G2 for 
2DG). We indicate the location uncertainties for these sources as 
circles in Fig.O We also show the locations of the powerful star 
cluster R 136 and of the two Crab-like pulsars PSR J0540-6919 
and PSR J0537-6910. From this figure it becomes clear that 
SI in the PS model does not coincide with any of these ob- 
jects. We confirm this observation by moving SI to the loca- 
tions of either PSR J0540-6919, PSR J0537-6910 or R 136, 
which leads to a significant TS reduction of 21.0 (4.2cr) for 
the position of PSR J0540-6919, 28.3 (4.9cr) for the position 



of PSR J0537-6910, and 31.1 (5.2o-) for the position of R 136. 
Thus, none of the 3 sources can explain the observed gamma-ray 
emission from 30 Dor alone. We come to the same conclusion 
when we replace G2 in the 2DG model by a point source at any 
of the 3 locations. 

To search for evidence of point source emission from 
PSR J0540-6919, PSR J0537-6910 or R 136 on top of the dif- 
fuse emission, we added point sources at the respective posi- 
tions to the 2DG model. Adding any of the sources led to only 
marginal TS improvements of 0.5 (0.3cr) for PSR J0540-6919, 
(Oct) for PSR J0537-6910, and 0.1 (O.lcr) for R 136. Adding 
all three sources together marginally improved the TS by 0.5 
(0.003cr). Thus, we found no evidence for point source emis- 
sion from any of the 3 sources on top of the diff'use emission 
that is modelled by G2. When we replaced G2 by a combination 
of 3 point sources located at the positions of PSR J0540-6919, 
PSR J0537-6910 and R 136 we obtained a TS of 1 1 13.4, which 
still is smaller by ATS = 9.2 than the TS obtained for the 2DG 
model. Obviously, the gamma-ray emission from 30 Dor is bet- 
ter described by a single 2D Gaussian shaped extended source 
than by a combination of point sources situated at the locations 
of PSR J0540-6919, PSR J0537-6910 and R 136. 

2.3.3. Gas templates 

To test specific hypotheses about the spatial distribution of the 
gamma-ray intensity we compared our data to spatial templates 
that trace the interstellar matter distribution in the LMC. The rea- 
soning behind this comparison is that gamma-rays are expected 
to primarily arise from interactions b etween cosmic rays and the 
interstellar medium (lOzel & Fichtellll 988), and we thus may ob- 
tain a measurement of the cosmic-ray density variations from the 
gamma-ray to gas ratio. 

Most of the gas in the LMC is found in form of neutral 
atomic hydrogen and helium (IStaveley-Smith et aL 1 12003.) while 
5% - 10% in mass is in form of molecular clouds dFukui et al.l 
2008) an d aboii t 1% in mass is in form of ionized hydrogen 
(Kennicutt et al. I [1995). For the atom ic hydrogen component 
(Hi) we used the map of Ki m et al.l (120051) that is based on 
a combination of aperture synthesis and multibeam data from 
ATCA and Parkes observations. We converted this map into hy- 
drogen column densities assuming optically thin emission or 
by applying an optical depth correction under the assumption 
of uniform spin temperature of either 40 K, 6 K or 125 K 
arx-Zimmer et al ] |2000h . [Bernard et al.l (l2008h suRRested the 
presence of an additional dark gas component in the LMC which 
correlates well with the distribution of the observed H i. Thus, 
if this dark gas component indeed exists it should be morpho- 
logically well described by the Hi templates. For the molecular 
hydrogen component (H2) we used th e CO map of the LMC ob- 
tained with the NANTEN telescope (iFukui et al. l 2008) which 
we converted into molecular hydrogen column densities using 
a constant CO intensity to H2 column d ensity conversion f actor 
Xco = 7 X 10^0 H2 cm-2 / (K km s'^) (iFukui et al.ll2008l) . For 
the ionized hydrogen component we used the full-sky composite 
Ha map of Finkbeiner ( 2003) which we converted into ionized 
hydrogen column densities using N(ll 11) = 1 .37 x IO^^/h^ which 
assumes a gas temperature of 8000 K and an ele ctron density of 
He = 1 cm"^ (cf. Eq. 6 of [Bernard et al .1120081) . In addition to 
the individual templates for the interstellar gas phases we also 
fitted a template of the total gas column density NiB) to the data 
that we derived by adding together the column densities of H i, 
H2, and Hii. Furthermore, we also fitted the linear combination 
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Table 4. Comparison of maximum likelihood model fitting re- 
sults (see text for a description of the models). 



Model 


TS 


Significance 


P 


PS 


1089.3 


31.5 


20 


2DG 


1122.6 


32.7 


10 


A^(Hl) (optically thin) 


771.8 


27.6 


2 


A^(Hi) (D+L components) 


792.9 


27.8 


4 


A^(Hl) {T, = 125 K) 


794.9 


28.1 


2 


A^(Hl) (r, = 60 K) 


834.6 


28.8 


2 


A^(Hl) (r, =40 K) 


835.3 


28.8 


2 


A^(H2) 


824.3 


28.6 


2 


A^(Hll) 


1110.1 


33.2 


2 


7V(H) 


806.9 


28.3 


2 


A^(Hl) + A^(H2) + A^(Hll) 


1110.7 


32.9 
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Note to the table: Columns are (1) the name of the model, (2) the TS of 
the fit, (3) the detection significance, and (4) the number of free param- 
eters. 



A^(Hi) + A^(H2) + A^(Hii) to the data with independent scaling 
factors and power law indices for each component. 

The TS values and detection significances that are obtained 
by maximum likelihood fitting are summarized in Tabled where 
we also give the corresponding results for the geometrical mod- 
els PS and 2DG for reference. Fitting the optically thin H i and 
H2 maps resulted in TS values that are considerably worse than 
those obtained for the geometrical models. Applying the optical 
depth correction to the H i data slightly improved the fit with a 
maximum TS of ~ 835 that is reached for spin temperatures of 
40 K - 60 K, yet still, this TS is considerably worse than that 
of the geometrical models. Apparently, the H i and H2 maps pro- 
vide rather poor fits to the data, indicating that the distribution of 
gamma rays does not follow the distribution of neutral hydrogen 
in t he LMC. 

iLuks & Rohlfsl (Il992h proposed that the neutral hydrogen 
in the LMC is confined into two separate structural features 
(dubbed 'D' and 'L' components), and we tested whether 
gamma-ray emission is possibly only correlated to one of these 
Hi components. We did this by separating the Hi data of 
iKim et all (l2005l) into two maps covering the heliocentric ra- 
dial velocity intervals 190 < Vrad < 270 km s"^ (L) and 270 < 
Vrad < 386 km s"^ (D) that roughly separate the two components 
(see Fig. 4 of iLuks & Rohlfsl 1 1992t) . Fitting both maps simulta- 
neously with independent scaling factors and power law indices 
to our data only marginally improved the fit upon the A/^(H i) gas 
map (ATS = 2L1). We therefore conclude that the gamma-ray 
emission does also not follow t he distribution o£either the D or 
the L component identified by Luks & Rohlfsl (Il992h . 

The H II map, on the other hand, gave a TS that is very close 
to that of the PS and 2DG models. Fitting the total hydrogen 
column densities, i.e. model A/^(H), gave a rather low TS, close 
to that obtained for the H i and H2 templates alone. This result 
is explained by the preponderance of atomic hydrogen in the 
LMC which means that the maps of total hydrogen and H i are 
very similar. Fitting a linear combination of the three gas phases, 
i.e. model N(R i) -h A/^(H2) + A/^(H 11), gave basically the same TS 
value that is obtained for the Hii template alone. In this case, 
Hii is in fact the only significant model component in the fit 
while the maximum likelihood fluxes for the other components 
are negligible. 

The H II map thus provided the best fit among all of the gas 
maps to the LAT data. Figure |4] shows the model counts map 



Table 5. Comparison of maximum likelihood model fitting re- 
sults after inclusion of the source G2 in the background model. 
Columns are identical to those in Table IH 



Model 


TS' 


Significance 


P 


A^'(Hi) (optically thin) 


173.1 


12.9 


2 


N'OAl) (D+L components) 


194.1 


13.7 


2 


A^'(Hi) {T, = 125K) 


175.4 


13.0 


2 


N'{m) (T, = 60 K) 


180.0 


13.2 


2 


N'iHl) (r, = 40 K) 


178.4 


13.1 


2 




197.0 


13.9 


2 


TV' (Mil) 


284.9 


16.7 


2 


N\R) 


183.7 


13.3 


2 


7V'(Hl) + 7V'(H2) + A^'(Hii) 


284.9 


16.1 
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for the H 11 template and the corresponding residual counts map, 
which apparently is similar to the residual maps that we obtain 
for the geometrical models. A search for point- like or extended 
emission on top of the H 11 template did not reveal any signifi- 
cant additional signal (ATS < 25). The H 11 map is characterised 
by a strong emission peak near 30 Dor which is attributed to the 
intense ionizing radiation of the massive stars in this highly ac- 
tive region. To test whether the fits are mainly driven by this 
peak of ionizing flux, we repeated the template map analysis 
by adding the source G2 which describes the emission from 
30 Dor in the 2DG model as an additional component to our 
background model. The flux and spectral index of this compo- 
nent were kept as free parameters. At the same time, we modified 
the gas maps by setting all pixels within a circular region of 0.5° 
around (0^12000, ^12000) = (05^38^42^ -69°06'03") to zero which 
removed any peak associated to 30 Dor from the templates (in 
particular, the bright peak of ionizing flux in 30 Dor is now re- 
moved). In that way, we considerably reduced the impact of the 
gamma-ray emission from 30 Dor on the fit of the gas templates. 

We summarise the results of this analysis in Table \5\ where 
the prime in the model names indicates that the 30 Dor emission 
has been removed from the gas templates. The TS values are now 
considerably reduced since Xo now includes the G2 source that 
accounts for gamma-ray emission from 30 Dor. Consequently, 
TS now measures the significance of gamma-ray emission that 
exists in addition to that seen towards 30 Dor in the LMC. To 
recognise this change in the definition we label the correspond- 
ing column in Table [5]by TS'0 Globally, however, we found the 
same trend as before: the Hii template provided the best fit to 
the data while the H i and H2 templates resulted in significantly 
smaller TS' values. We thus conclude that the morphology of 
the gamma-ray emission from the LMC is not well explained 
by the distribution of neutral gas in that galaxy, irrespectively 
of whether the emission from 30 Dor is taken into account. A 
template based on the distribution of ionized hydrogen as mea- 
sured by the Ha emission describes the spatial distribution of 
the gamma-ray emission in the LMC considerably better than 
the distribution of neutral gas. 

2.4. Variability 

As we noted earlier (cf. section 12.21) , our data indicated possi- 
ble flaring activity within the LMC boundaries during month 4 



^ Fitting the G2 source as the only component to the data resulted in 
TS= 826.5. We can thus formally convert TS' into TS using TS=TS' + 
826.5. 
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Fig. 5. Monthly lightcurve for gamma-ray emission from the 
LMC obtained using the H ii gas map. The 4th month of data that 
was excluded from our analysis is indicated by an open symbol 
and grey colour. The solid line indicates the average flux over 
the full 1 1 months of data that is obtained by assuming a power 
law spectral shape for the emission. 



(MJD 54777.8 - 54808.2) of the observation interval, and we 
excluded the data of this month from our analysis. Now with 
reliable spatial templates for the LMC gamma-ray emission at 
hand (the 2DG model and the H ii gas map), we reconsidered the 
time variability of the gamma-ray emission in the LMC area. 

First, we considered the time-variability of the integrated 
gamma-ray emission from the LMC. For this purpose we fitted 
the 2DG model and the H ii gas map to the data on a monthly, 
2 weeks and weekly basis using a power law to describe the 
source spectrum. We fixed the power law spectral index to the 
average value that we obtain from fitting the 2DG model or the 
Hii gas map to all data excluding month 4. Figure [5] shows the 
monthly lightcurve that we obtain using the H ii gas map. Fitting 
the lightcurve with a constant flux level (solid line) resulted in 

= 25 A for 10 degrees of freedom, which corresponds to a 
probability of 0.2% for the flux being constant. Obviously, the 
flux measured during month 4 (MJD 54777.8 - 54808.2) from 
the LMC is significantly larger than the flux measured for the re- 
maining period. Similar results were obtained when we perform 
the variability analysis on time scales of 2 weeks (x^ = 38.9 for 
23 degrees of freedom, L0% probability for constant flux) and 
1 week (x^ = 69.9 for 47 degrees of freedom, 0.8% probability 
for constant flux). Using the 2DG models instead of the Hii gas 
map gave comparable results. 

As next step we tried to localise any flaring source (and in 
particular the one that occurred during month 4) by fitting a 
point source with free position, flux and spectral power law in- 
dex on top of the 2DG model and the Hii gas map to each of 
the 11 monthly datasets. We fixed the flux and spectral index 
of the 2DG model and H ii gas map to the values that we ob- 
tained from the fit to all data excluding month 4. This analysis 
revealed that month 4 is the only period where a point source 
is detected significantly (TS > 25) on top of the LMC mod- 
els within the galaxy's boundaries. Using the Hii map as tem- 
plate for the LMC, we localised the source at (aj2ooo, ^j20oo) = 
(05M6.4"^, -69°0r) with a 95% containment radius of 29' and a 
detection significance of 4.6cr (TS = 30.1), while with the 2DG 
model we found (aj2ooo, ^j2ooo) = (05M6.2"^, -69°0r) with a 
95% containment radius of 36', and a detection significance of 
4.5o- (TS = 29.3). Obviously, the choice of the LMC model had 



a negligible impact on the maximum likelihood location of the 
source. 

We show the larger localisation error circle obtained with the 
2DG model in Fig. [3l which illustrates that the flaring source 
is located near 30 Dor and that it is close to the maximum 
gamma-ray intensity observed from the LMC. Independent on 
the LMC model, we found an integrated > 100 MeV photon flux 
of (21 ±8 (stat)±4 (sys))x 10"^ ph cm"^ s"^ for the flaring source 
during month 4 and obtained a maximum likelihood power law 
spectral index of F = 2.8 ± 0.3 (stat) ±0.1 (sys). The spec- 
tral index is rather steep and compatible with the softest indices 
that are typically found for flat- spectrum radio quasars (FSRQs) 
( Abdo et al. 2009b). FSRQs tend to show also the l argest vari- 
abilities among all source classes detected by LAT (lAbdo et al.l 
i2009aiib) which makes them good candidates for flar ing sources. 
Using t he source association procedure described by I Abdo et al.l 
(l2009ah . we searched for plausible counterparts of the flaring 
source. We did not find a plausible counterpart in any of our 
standard blazar or AGN catalogues. The only object that ful- 
filled our association criterion was RX J0546. 8-6851 which is 
listed in the Magellanic Cloud high-mass X-ray binary catalogue 
of iLiu et al.l (|2005). However, the nature of t his source that has 
been discovered by the Einstein Observatory (iWang et al.ll 19911) 
is not yet established. So far, no compani on star has been identi- 
fied for RX J0546.8-6851, and following [Kahabka et al.l (l200lh 
the source could also be an active galactic nucleus. It is thus con- 
ceivable that the flaring source detected by LAT during month 4 
is not physically associated with the LMC, but is rather a yet 
unidentified background blazar that by chance is located near 
the line of sight toward 30 Dor. 

2.5. Blazars in the field of the LMC 



As we noted earlier (cf. section 12.21) , we expect to find 1-2 
gamma-ray blazars within the boundaries of the LMC in our 
dataset. Identification of these blazars would be possible if they 
were flaring during the observation period (see the previous sec- 
tion), yet a substantial fraction of the blazars seen by the LAT 
do not show any significant time variability ( Abdo et al.ll2009bh . 
Blazars are thus difficult to recognize if they spatially concur 
with gamma-ray emission from the LMC. Nevertheless, we can 
try to spot blazar candidates by searching for point source emis- 
sion that spatially coincides with the locations of blazars. 

For this purpose we searched for counterparts for any of the 
point sources in t he PS model. Using our standard source asso- 
ciation procedure (lAbdo et al.ll2009al) we did not find any blazar 
candidate coincident with any of the sources of the PS model. 

As next step, we simultaneously fitt ed point soii r ces at 
the positions of all CRATES blazars (iHealev et al.l l2007h 
that are found within a radius of 5° around the centre of 
the LMC, which we assumed here to be (0^12000, ^j200o) = 
(05M7^36^ -69°0r48"), in addition to the 2DG model and the 
H II template to the data. The fluxes and the spectral power law 
indices of all sources were left as free parameters of the fit. 
We found a total of 36 CRATES sources within the selection 
region of which one (CRATES J060106-703606) was already 
part of our background model (see Table [T]). This left us with 
35 CRATES blazars in the field. As expected from our previ- 
ous searches for significant point sources above the 2DG model 
and the Hii template (cf. sections [2.3.21 and |2 . 3 . 3K none of these 
blazars was significantly detected (TS > 25). However, two 
sources had a formal detection significance above 2cr^ CRATES 



which corresponds to TS ^ 6 for p -2. 
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Fig. 6. Gaussian kernel (cr = 0.2°) smoothed residual map for 
the 2DG model over which we overplot the CRATES blazars in 
the field of the LMC as stars. The boxes mark those blazars for 
which TS > 6. 



105 1636-723707 (TS = 16.7), and CRATES J053524-673430 
(TS = 12.4). The locations of these blazars are overplotted to- 
gether with those of all other CRATES sources in the LMC field 
over the residual map for the 2DG model in Fig. (6] 

CRATES 105 1636-723707 app eared to be point like and had 
already been associated in section 12.3.21 as plausible counter- 
part of the 3rd 2D Gaussian in the 2DG model analysis. Figure 
[6] clearly shows that this source is associated with a peak in 
the observed counts maps. CRATES J053524-673430, in con- 
trast, is not associated with any emission peak but is rather sit- 
uated in the centre of an extended emission feature. Fitting a 
2D Gaussian shaped emission profile at the location of CRATES 
J053524-673430 suggested indeed that the emission is extended 
(cr = 26' ± 10'), making the blazar an unlikely counterpart of the 
gamma-ray emission in that area. CRATES J05 1636-723707 
seems to be thus the only plausible blazar candidate in the field 
of the LMC. 

We note, however, that although many of the LAT high- 
latitude sources have counterparts in the CRATES catalogue, 
there is still a non-negligible population of unassociated high- 
latitude sources in the LAT source list that may consist of yet 
unknown (or unidentified) blazars. Some of the residuals seen in 
Fig. [6] could indeed arise from such blazars, yet their emission 
does not appear to significantly contribute to the overall emis- 
sion from the LMC region. 

2.6. Pulsar Hghtcurves 

As we noted earlier (cf. section 12.31 and Fig. [3]), the peak of 
gamma-ray emission that we see toward 30 Dor is spatially close 
to the locations of the powerful young pulsars PSR J0540-6919 
and PSR J0537-6910, which both are potential sources of high- 
energy gamma rays. By searching for their pulsations in our data 
we can assess if the pulsars indeed contribute to the observed 
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Fig. 8. Spectrum of the LMC emission obtained by fitting spa- 
tial model maps in 6 logarithmically spaced energy bins covering 
the energy range 200 MeV - 20 GeV. Black dots show the total 
spectrum obtained by fitting the H n gas map, red squares and 
blue triangles show the spectra for the LMC disk and for 30 Dor, 
respectively, obtained using the 2DG model. Error bars include 
statistical and systematic uncertainties. The lines show the re- 
sults obtained by fitting exponentially cut off' power law spectral 
models to the data using maximum likelihood fits. 



emission. We used for this search ephemerides that have been 
obtained from concurrent RXTE observations (Marshall, private 
communication) and that cover the time intervals MJD 54651- 
55015 for PSR J0540-6919 and MJD 54710-54714, 54751- 
54771, and 54885-54897 for PSR J0537-6910. Only data within 
these time intervals have been used for the analysis. Gamma-ray 
photon arrival times were referred to the solar-system barycentre 
and pulse phases w ere assigned using t he standard pulsar timing 
software TEMP02 (iHobbs et al.ll2006h . 

Figure [7] shows the Hghtcurves that we obtained by select- 
ing photons within a radius of 0.5° around the nominal pul- 
sar locations for the energy range 200 MeV - 20 GeV. For 
PSR J0540-6919 our data may indicate the presence of pulsa- 
tions, but the H-test statistic amounts to only 10.1 which cor- 
responds to a detection significance of 2 Act for the pulsations. 
For PSR J0537-6910 the lightcurve is featureless and an H- 
test gives 3.2 corresponding to a detection significance of l.lcr. 
Thus, neither pulsar is detected significantly in our da ta which 
is in line with our previous finding (cf . section 12.3.2b that the 
gamma-ray emission from 30 Dor cannot be satisfactorily de- 
scribed by point source emission from the pulsars. 

2.7. Spectrum 

So far, all analysis has been done by assuming that the spectrum 
of the LMC emission components is well described by power 
laws. To determine the spectrum of the gamma-ray emission 
from the LMC independently of any assumption on the spec- 
tral shape, we fitted our data in 6 logarithmically spaced energy 
bins covering the energy range 200 MeV - 20 GeV. We obtained 
the total spectrum of the LMC by fitting the H n template to the 
data. We also obtained separate spectra for the LMC disk and for 
30 Dor by fitting the 2DG model to the data. Here, Gl is taken to 
represent the LMC disk emission while G2 is taken to represent 
the emission from 30 Dor. The results are shown in Fig.O The 
spectra are relatively fiat below a few GeV (F ~ 2) and show 
evidence for a break or cut off' above that energy. 
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Parameter 


Total 


Disk 


30 Doradus 


Photon flux (>100 MeV) 


26.3 + 2.0 + 4.2 


18.4 + 2.3 + 3.0 


6.9 + 1.7 + 0.7 


Energy flux (>100 MeV) 


16.0 + 0.6 + 1.8 


11.1 +0.7 + 1.5 


4.6 ± 0.5 + 0.5 


Photon flux (>200 MeV) 


12.2 + 0.4 + 0.7 


8.3 + 0.6 + 0.7 


3.3 + 0.4 + 0.3 


Energy flux (>200 MeV) 


12.9 + 0.4+1.8 


8.9 + 0.5 + 1.6 


3.8 + 0.3 + 0.5 


r 


2.04 + 0.09 + 0.07 


2.10 + 0.13 + 0.10 


1.96 + 0.25 + 0.02 


(GeV) 


7 1 +3.5+6.6 
^•-^-2.1-2.4 


10 3 + 12-4+00 
-^^•-^ -4.4- 10.3 


(1 Q+lO.5+1.3 
'^•^-3.5-1.0 


TS with respect to power law 


11.3 


3.0 


3.6 


Significance of cut off (cr) 


3.4 


1.7 


1.9 



Note to the table: Photon fluxes are given in units of 10 ^ ph cm ^ s \ energy fluxes are given in units of 10 erg cm ^ s ^ Errors are statistical 
and systematic. 



To determine the integrated gamma-ray flux of the spectra 
we fitted exponentially cut off' power law spectral models of the 
form N{E) = k (E/Eq)'^ Qxp(-E/Ec) to the data. We made these 
ffts by means of a binned maximum likelihood analysis over the 
energy range 200 MeV - 20 GeV. This analysis is more reli- 
able than fftting the spectra of Fig. [8] directly since it accounts 
for the Poissonian statistic of the data. The resulting spectral 
shapes are shown as lines in Fig.[8]and spectral fft parameters are 
given in Tabled Integrated ffuxes are derived for energies above 
100 MeV and 200 MeV by integrating the spectral model numer- 
ically from the lower boundary up to 500 GeV. The > 100 MeV 
ffuxes are mainly given for comparison with the former measure- 
ments done by EGRET. Since we restricted our analysis to pho- 
tons > 200 MeV the extrapolation down to 100 MeV introduces 
an additional systematic uncertainty that for the photon ffux can 
reach up to 30% (cf. Table [S]). We also determined the signif- 
icance of the spectral cut off by computing the TS difference 
between fitting a simple power law and fitting an exponentially 
cut off power law for each component. 

3. Discussion 

3.1. Comparison to EGRET 

ISreekumar et al.l (Il992k reported the first detection of the LMC 
in high-energy gamma rays based on 4 weeks of data collected 
with the EGRET telescope aboard CGRO. They obtained an in- 
tegrated > 100 MeV ffux of (1.9 ± 0.4) x 10"^ ph cm'^ s"^ for 
the LMC from their data which is in reasonable agreement with 
our value of (2.6 ± 0.2 (stat) ± 0.4 (sys)) x 10"^ ph cm'^ s"^ that 

we obtained from spectral mode l fftting. 

The 3rd EGRET catalogue (iHartman et al.lll999! ) gives po- 
sitions, ffuxes and power law spectral indices for point sources 
detected during 4 years of EGRET observations and locates the 
LMC emission at (aj2ooo, ^J20oo) = (05^33"^, -69° 160 with a 
95% containment radius of 32\ close to our maximum likeli- 
hood location for S 1 in the PS model or G2 in the 2DG model, 
and consistent with the location of 30 Dor. The > 100 MeV pho- 
ton ffux in the 3rd EGRET catalogue amounts to (1 .4±0.2)x 10"^ 
ph cm"^ s"^ which is signiffcantly lower than our value. The 
assumption of a single point source for the EGRET catalogue 
analysis may have led to an underestimation of the ffux from 
the extended LMC galaxy. On the other hand, our > 100 MeV 
ffux is based on power-law extrapolation from our measurement 
> 200 MeV which could lead to overestimation at lower en- 
ergies. The 3rd EGRET catalogue quotes a power law spectral 
index of F = 2.2 ± 0.2 for the LMC which is consistent with our 
value. In particular, fftting the total LMC emission with a simple 
power law instead of an exponentially cut off power law we ob- 



tain F = 2.26 ± 0.04 (stat) ± 0.10 (sys), consistent with the value 
obtained from EGRET data. 

Thus, globally, our results are in line with those found earlier 
by EGRET, yet we benefft from a much better photon statistic^, 
an improved angular resolution and spectral coverage up to at 
least ~ 10 GeV. 

3.2. On the origin of gamma rays from 30 Doradus 

Although it was already obvious from EGRET observations that 
the LMC was an extended source, only the Fermi/LAI data al- 
low now for the ffrst time to clearly resolve the LMC and to 
localise the emission maximum near the 30 Dor massive star 
forming region. Even if we cannot establish either of the pul- 
sars PSR J0540-6919 or PSR J0537-6910 as th e sour ce of the 
gamma-ray emission toward 30 Dor (cf. sections [2.3.21 and [231) , 
the indication of possible pulsations from PSR J0540-6919 is 
a challenging possibility. If further observations of the LMC by 
Fermi/LAI conffrm the detection of pulsations, this would be 
the ffrst time that a pulsar outside our own Galaxy were de- 
tected in gamma rays. We are aware that the H.E.S.S. collab- 
oration has announcecfl the detection of the pulsar wind nebula 
of PSR 105 37 -69 10 at TeV energies, yet so far we do not ffnd 
evidence for gamma-ray emission from that source in our data. 

As the gamma-ray emission from 30 Dor is not dominated 
by pulsars, cosmic-ray interactions with the ambient interstel- 
lar medium and radiation ffelds are plausibly the origins of the 
observed emission. In this case, the small size of the emis- 
sion region puts stringent constraints on the diffusion length 
of GeV cosmic-ray protons. For source G2 we ffnd an extent 
of (T = IT ± 4' and a 95% conffdence level upper limit of 
cr = 19\ which at the distance of the LMC corresponds to 
linear dimensions of 170 pc and 270 pc, respectively. This is 
extremely small compared to expectations. Gamma-ray obser- 
vations of our own Galaxy suggest in fact a Galactocentric 
cosmic-ray emissivity gradient that is much smaller than that 
of the distribution of supernova remnants or pulsars, which has 
been taken as evidence for large-s ca le cosmic-ray diffusion in 
our G alaxy (Bloemen etaD [I984allbl [19861: [Strong et al.[ [1988[: 
[Stro ng & Mattox 1996). In addition, the low average gas density 
that is encountered by cosmic rays as determined from cosmic- 
ray isotope ratios indicates that the particles spend a consider- 



^ The total number of counts we attribute to the LMC amounts to 
~ 1550 photons while the 3rd EGRET catalogue quotes a total of 192 
counts obtained during 4 years of observations. 

^ The detection has been announced at the 3 International Cosmic- 
Ray Conference in Lodz 2009, however no details about source flux and 
spectrum are yet publicly available. 
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able fraction of their lifetimes outside the plane of the Galaxy, 
suggesting a cosmic-ray halo that extends to scale heights of 
4-6 kpc above t he galactic plan e (Moskalenko et al. 2004). On 
the other hand, iFichtel & KniffenI (119841) have shown that the 
gamma-ray data are also c ompatible with a tight coupling of cos- 
mic rays to the matter, and lStrong et alJ (l2004bl) suggest that ne- 
glecting a possible strong Galactocentric gradient in the CO-to- 
H2 conversion factor has led in the past to an underestimation of 
the Galactic cosmic-ray gradient. Thus, observations may also 
be compatible with a short cosmic-ray diffusion length in our 
own Galaxy. 

Galactic studies are however hampered by our position 
within the Galactic plane, which does not allow for an unam- 
biguous correlation of the gamma-ray emissivity with Galactic 
star forming regions, making estimates of cosmic-ray diffusion 
rather uncertain. Our observations of the LMC provide for the 
first time a much clearer picture about the link between possible 
cosmic-ray acceleration sites and gamma-ray emission, and sug- 
gest that both are tightly connected. 30 Dor would thus be the 
most powerful cosmic-ray acceleration site in the LMC, which 
would be no surprise as it is one of the mostactive starburst re- 
gions in the entire Local Group of galaxies dPineda et al.ll2009i) . 

3.3. Gamma rays from cosmic-ray interactions 

The gamma-ray emission from cosmic-ray interactions with the 
interstellar medium and radiation field is expected to be in- 
trinsically diffuse in nature. Indeed, our analysis suggests that 
the gamma-ray emission from the LMC does not originate in 
a small number of individual point sources. By adding succes- 
sively point sources to our model (cf. section [TT2b we formally 
detected 5 significant point sources (model PS), but we found 
a model with less components that results in a higher TS value 
when we allow for the source components to be extended (model 
2DG). It is thus more likely that the LMC emission is indeed dif- 
fuse in nature, or alternatively, composed of a large number of 
unresolved and faint sources that can not be detected individu- 
ally by Fermi/LAI. 

From here on we thus make the assumption that the gamma- 
ray emission from the LMC indeed originates from cosmic-ray 
interactions with the interstellar gas and the stellar radiation 
field, and we discuss the implications of this assumption in light 
of our observations. In particular we test the 2 extreme hypothe- 
ses that (1) all gamma-ray emission from the LMC is attributed 
to diffuse emission from cosmic-ray interactions, and (2) only 
the emission from component Gl in the 2DG model arises from 
cosmic-ray interactions, while the gamma rays from 30 Dor (or 
component G2 in the 2DG model) originate from other sources. 
Hypothesis 2 is mainly motivated by the particular nature of 
30 Dor, and we want to examine how our conclusions are altered 
when we exclude this source from our considerations. Note that 
in this case, only 70% of the integrated LMC flux is attributed to 
cosmic-ray interactions while for hypothesis 1 we assume that 
100% of the LMC flux is due to cosmic rays. 

Using these two hypotheses (which we call from here on 
'Hi and 7^2) we computed the average integrated > 100 MeV 
gamma-ray emissivity per hydrogen atom of the LMC using 



where is the integrated > 100 MeV photon flux from the 
LMC under "Hi or 7^2, and J A^h is the spatially inte- 
grated hydrogen column density of the LMC. From the H i map 



Table 7. Results of the gamma-ray emissivity analysis of the 
LMC for hypotheses "Hi and 'H2. 



Parameter 


•Hi 


^2 




25.0+1.7 + 4.1 


15.2 + 1.9 + 3.4 


qy 


6.9 + 2.4 + 2.6 


4.2 + 1.5 + 1.7 


rc 


0.31+0.01+0.02 +^13 


0.21 + 0.01 + 0.02 +0-^9 



Note to the table: Given are the integrated > 100 MeV photon flux 
in units of 10~^ ph cm~^ s~\ the average > 100 MeV gamma-ray 
emissivity qy in units of 10~^^ ph s~^ sr~^ H-atom~\ and the cosmic-ray 
enhancement factor r^. Quoted errors are statistic and systematic. For 
rc we quote in addition the (asymmetric) systematic error related to the 
uncertainty in the total LMC gas mass. 

of iKim etal.1 (l2005l) one obtains / Vr dQ = 2.39 x lO^^ H- 
atom cm"^ sr for atomic hydrogen, assuming optically thin emis- 
sion. Including molecular hydrogen increases the value by 5% - 
10% and applying optical depth corrections for 125 K and 60 K 
further increases the integral by 12% and 42%, respectively. 
Evidence for some hidden hydrogen (or the presence of substan- 
tial optically thick emission) comes from an analysis of infrared 
data obtained with Spitzer, suggesting an additional component 
that follows the distribution of the H i seen in the radio surveys, 
and that doubles the total Ras mass (iBernard et al.ll2QQ8i) . Thus, 
J Vh dQ could be as large as 4.8 x 10^^ H-atom cm~^ sr and 

we adopted J Vr dQ = (3.6 ± 1.2) x 10^^ H-atom cm"^ sr as a 
mean value and included the uncertainty of the total hydrogen 
mass of the LMC in our computation. The resulting emissivities 
are given in Table [71 

Depending on the hypothesis we make, the average gamma- 
ray emissivity of the LMC is between ~ 2 to ~ 4 times smaller 
than the > 100 MeV emissivity of qy = (1.63 ± 0.05) x lO'^^ 
ph s"^ sr"^ H-atom"^ that has been determined by Fermi/LAI 
for the local interstellar medium of our own Galaxy (lAbdo et al.l 
l2009c l). We illustrate this difference in Fig. |9j where we com- 
pare the differential gam ma-ray emis sivity of the local interstel- 
lar medium (cf. Fig. 5 of 'A bdoet all 2009c) to that of the LMC. 
Following Eq. ([T]), we derived the differential gamma-ray emis- 
sivity by dividing our spectra by J Vh dD = (3.6 ± 1 .2) x 10^^ H- 
atom cm"^ sr. For *Hi we used the spectrum derived using the 
Hii gas map for the total emission from the LMC (black dots), 
while for *K2 we used the LMC disk spectrum obtained using 
component Gl of model 2DG (red dots) 

We compared the differential gamma-ray emissivities to a 
one-zone model of cosmic-ray interactions with the interstel- 
lar medium that takes into account decay following proton- 
proton interactions, Bremsstrahlung from cosmic-ray electrons 
and inverse Compton scattering of cosmic-ray electrons on LMC 
optical and infrared photons and cosmic microwave background 
photons. We calculated the tt^ production by proton-proto n in- 
teractions following the prescription of K amae et al.l (l2006l) . We 
used the cosmic-ray proton, e lectron and positron spectra pre- 
sented in lAbdo et all (l2009d) for the local Galactic environ- 
ment. The scaling factor (r^) of our model with respect to the 
data, which we refer to as the cosmic-ray enhancement factor, 
is thus a direct measure of the average LMC cosmic-ray den- 
sity with respect to that in the vicinity of the Earth. The tt^ 
and Bremsstrahlung emissivities were calculated by assuming 
the average LMC metallicity to be Z = 0.4 Z© (WesterluncJ 
1 19971) . We applied a corresponding nucl ear enhance ment fac- 
tor of 6m = 1.75 to the tt^ emissivities (lMorill2QQ9l) and took 
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Fig. 9. Differential average gamma-ray emissivity for *Hi (black 
dots) and 7^2 (red dots) compare d to that of the Galac tic local in- 
terstellar medium determined by Abdo et aP (l2009c) (grey data 
points). The error bars for "Hi and include statistical and sys- 
tematic errors and the uncertainties in the total gas mass. The 
solid line shows the predicted gamma-ray emissivity computed 
in the framework of a one-zone model for (see text). The 
other lines show the contributions of tt^ -decay (long dashed), 
Bremsstrahlung (short dashed), and inverse Compton emission 
(dotted). 

the total hydrogen mass of the LMC to be 7.2 x 10^ Mo0 The 
inverse Compton co mp onen t w as calcula ted using the method 
described by iBlumenthal & Gouldl (|197Q|) using the optica l and 
infrared interstellar radiation fields from lPorter et al.l (12008 ) that 
we rescaled according to th e stellar luminosity density and the 
observed infrared emission (Berna rd et al.l l2008V respectively. 

Fitting this spectral model to our data using a binned max- 
imum likelihood analysis gave an average cosmic-ray enhance- 
ment factor of rc = 0.31 ± 0.01 for and of = 0.21 ± 0.01 
for 7^2 (cf. Table O. Systematic errors due to uncertainties in the 
eff'ective area of the instrument amount to ±0.02. An additional 
systematic error of -23% to -1-42% comes from the uncertainty 
in the total gas mass of the LMC, which largely dominates the 
statisti cal and sy s temat ic measurements errors. 

Fich tel et aP dlQQlh estimated the cosmic-ray enhancement 
factor to be Tc ~ 1 by comparing the synchrotron radiation of the 
LMC to that of the Milky Way (cf. their Fig. 3). We also com- 
puted the synchrotron emission in our model under the assump- 
tion that the cosmic-ray proton to electron ratio in the LMC and 
the electron spectrum are the same th an th ose measured in the 
local interstellar medium (lAbdo et al.ll2009ch . We further took 
the thickness of the emi tting region to be 2 kpc and assumed a 
magnetic field of 5 yuG (|Pohllll993h . Our predicted synchrotron 
flux is in good agreement with th e observ ed non-thermal flux 
quoted by Klein et al. (1989). Fich tel et al.1 (l991), however, as- 
sumed a highe r non-thermal flu x for the LMC compared to the 
observations of lKlein et al.l (Il989l) . which may explain why their 
Vc i s in excess of our r n odel. 

Pavlido u & Fieldsl (1200 ll) followed another approach and 
predicted = 0.14 from a comparison of the supernova rate 
in the LMC to that in our Galaxy, yet they considered also 
higher values due to the large uncertainties that are inherent 
to the estimates of supernova rates. Their estimate relies fur- 
ther on the assumption that the cosmic-ray containment in the 

^ For a distance of 50 kpc to the LMC Ja^h dQ = (3.6+1.2)xl0i9 H- 
atom cm~^ sr corresponds to a total H i mass of (7.2 + 2.4) x 10^ M©. 



LMC is similar to that in our own Galaxy, which is not nec- 
essarily fulfilled since cosmic rays may escape fairly rapidly 
into intergala ctic space because of the small size of the LMC 
(iGinzburg & P tuskin 1985). Our value, which is slightly above 
the estimate of Pavlidou & Fieldsl (1200 Ih . suggests however that 
cosmic-ray escape from the LMC is not more important than for 
our own Galaxy. This is in line with the relatively small proton 
diff'usion length that we inferred from the compactness of the 
30 Dor source (cf. section [T2b . 

3.4. The sites of cosmic-ray acceleration in the LMC 

To reveal the sites of cosmic-ray acceleration in the LMC we 
mapped the cosmic-ray density variations in the galaxy by com- 
puting the gamma-ray emissivity qy as function of position. We 
did this by generating a background subtracted counts map for 
and a background and G2 model subtracted counts map 
for *H2 that we divided by the N(tl) map after convolution of 
the latter with the LAT instrumental response function. We nor- 
maUsed A^(H) to a total LMC hydrogen mass of 7.2 x 10^ Mq 
that takes into account the possible presence of dark gas that 
is not seen in ra dio surveys of H i ( see section [33l) . We adap- 
tively smoothed (lEbeling et al.ll2006h the counts maps and used 
the resulting smoothing kernel distribution to smooth also the 
convolved A/^(H) map before the division to reveal significant 
structures at all possible scales, while suppressing the noise that 
arises from the limited photon counting statistics. 

The resulting emissivity maps are shown in Fig. \TU[ We 
superimpose on the images the interstellar gas distribution, as 
traced by A/^(H), convolved with the LAT instrumental response 
function, and also show the locations of potential particle accel- 
eration sites, such as pulsars, supernova remnants, Wolf-Rayet 
stars and supergiant shells. 

Figure [TOl reveals that the cosmic-ray density varies consid- 
erably over the disk of the LMC. The gamma-ray emissivity is 
highest in 30 Dor and the northern part of the galaxy, while the 
southern part and in particular the dense ridge of gas south of 
30 Dor seems basically devoid of cosmic rays. The large varia- 
tions in qy confirm our earlier findings (cf. section [2331) that the 
gamma-ray emission correlates little with the gas density in the 
LMC, and this is irrespective of whether we consider 30 Dor to 
be part of the emission (left panel versus right panel). The most 
striking variation is found along the gas ridge that runs over ~ 3° 
along aj2ooo ~ 05^40"^ which coincides with the most prominent 
region of ^^C O emissio n tracing giant molecular clouds in the 
LMC (Fu kui et al.lfT999 ). Roughly 2 0% of the total gas m ass in 
the LMC is confined into this ridge ( iLuks & Rohlfslll992h . and 
if the cosmic-ray density were uniform over the LMC (or at least 
over the ridge area) the entire ridge should be a source of high- 
energy gamma rays. This, however, is obviously not the case. 
Cosmic rays apparently do not penetrate into the southern part 
of the ridge, which is an additional argument in favour of a short 
GeV cosmic-ray proton diff'usion length in the LMC. 

Figure [TOl suggests further that the cosmic-ray density corre- 
lates with massive star forming tracers, and in particular Wolf- 
Rayet stars and supergiant shells. This correlation is particularly 
striking if we assume that the emission from 30 Dor is originat- 
ing in the same physical process as the rest of the LMC emission 
(*Ki). If we exclude 30 Dor from the comparison CH2) the corre- 
lation becomes less striking, yet still, most of the cosmic rays are 
found in regions that are rich in Wolf-Rayet stars and that are lo- 
cated near supergiant shells. This finding is corroborated by the 
good fit of the H 11 gas map (which is probably the most direct 
tracer of massive star forming regions within a galaxy), even 
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if the 30 Dor emission is discarded from the analysis (cf. sec- 
tion [233l) . 

Thus, the gamma-ray emissivity maps of the LMC support 
the idea that cosmic rays are accelerated in massive star forming 
regions as a result of the large amounts of kinetic energy that are 
input by the stellar winds and supernova explosions of massive 
stars into the interstellar medium. Our data reveal a relatively 
tight confinement of the gamma-ray emission to star forming 
regions, which suggests a relatively short diff'usion length for 
GeV protons. 

4. Conclusions 

Observations of the LMC by Fermi/LAI have for the first time 
provided a detailed map of high-energy gamma-ray emission 
from that galaxy. Our analysis revealed the massive star form- 
ing region 30 Doradus as bright source of gamma-ray emis- 
sion in the LMC in addition to fainter emission regions found 
in the northern part of the galaxy. So far, we could not identify 
any point source contributions to the emission from the 30 Dor 
region, and in particular we did not significantly detect pul- 
sations from the energetic Crab-like pulsars PSR J0540-6919 
and PSR J0537-6910, although the lightcurve of the former 
indicates variability at the 2Acr level. Cosmic-ray interactions 
with the interstellar medium and radiation field seem thus the 
most plausible origin for the observed gamma-ray emission from 
30 Dor. 

The gamma-ray emission from the LMC shows very little 
correlation with gas density. A much better correlation is seen 
between gamma-ray emission and massive star forming regions, 
as traced by the ionizing gas, Wolf-Rayet stars and supergiant 
shells, and we take this as evidence for cosmic-ray acceleration 
in these regions. The close confinement of gamma-ray emission 
to star forming regions implies a relatively small GeV cosmic- 
ray protons diff'usion length. 

Continuing observations of the LMC with FermifLAI in the 
upcoming years will provide the photon statistics to learn more 
about the origin of the gamma-ray emission from that galaxy. 
Better statistics will help in identifying more individual emission 
components and may help to separate true point sources from the 
more diff'use emission that we expect from cosmic-ray interac- 
tions. In particular, 1-2 year's of additional exposure should al- 
low the question of the contribution of PSR J0540-6919 to the 
gamma-ray emission from 30 Dor to be answered, and should 
also allow clearly establishing the pulsar wind nebulae contribu- 
tions at higher energies. 
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Fig. 4. Gaussian kernel (cr = 0.2°) smoothed model counts maps (left panels) and residual maps (right panels) of the region of 
interest (ROI). The model counts maps have been obtained by convolving the model sky maps that were scaled to the fitted flux 
with the LAT point spread function. The top row shows the PS model, the middle row shows the 2DG model and the bottom row 
shows the Hii template. The colour scaling and the dynamic range for the residual maps is identical to the right panel of Fig. [T] to 
allow for comparison. 
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LAT collaboration: Observations of the Large Magellanic Cloud with Fermi 

PSR J0540-691 9 PSR J0537-691 
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Fig. 7. Lightcurves obtained for PSR J0540-6919 (left) and PSR J0537-6910 (right) for the energy range 200 MeV - 20 GeV 
using an event selection radius of 0.5° around the nominal pulsar positions. The data are binned into 10 phase bins. The dotted 
lines indicate the sum of the expected contributions from extragalactic and Galactic diffuse emissions and residual instrumental 
background. 
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Fig. 10. Integrated > 100 MeV emissivity maps of the LMC in units of 10"^^ ph s"^ sr"^ H-atom"^ for *Hi (left panel) and 
(right panel). An adaptive smoothing with a signal-to-noise ratio of 5 has been applied to reduce statistical fluctuations. Note that 
the colour scales are not the same for the two images. As contours we show N(ll) column densities convolved with the LAT point 
spread function that are linearly spaced from 1 0% - 90% of the m aximum in steps of 10%. Symbols indicate the locations of pulsars 
(pluses) from ATNF catalogue version 1.36 ("Manchester et al.ll2Qb5i) , supernova remnants (diamonds) from Ro sa Williams web 
page http://www.astro.illinois.edu/projects/atlas/index.ht ml, Wolf-Raye t stars (stars) from the fourth catalogue of iBreysacher et al.l 
(119991) , and supergiant shells (circles) from lStaveley-Smith et al.l (l2QQ3l) . 



